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I. INTRODUCTION 



The hadron structure is believed to be described in terms of the fundamental theory of strong interactions, quantum 
chromodynamics (QCD), whose equations are notoriously difficult to solve. A successful approach in high-energy scat- 
tering is based on light-front quantization where hadrons are described by light-cone wave functions (LCWFs) [l| . The 
,^3 ' latter are expressed as an expansion of various quark (g), antiquark (q) and gluon (<?) Fock components. Schematically, 
Or a nucleon state is conceived as the following superposition 

Oh' 

CD ' \n) = V(3<r) mi) + i>{3q+ig)mm + tp(3q+qq)\qqqqq) + (i) 
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~ ~ '■ where in the light-cone gauge A + = the valence three-quark LCWF VY3g) involves six independent amplitudes 
| corresponding to different combinations of quark orbital angular momentum and helicity p| , and the Fock component 
^(3g+ig) with three quarks plus one gluon involves 126 independent amplitudes [|[. Adding a pair of sea quarks into 
00 ' the valence component to build the amplitude ^(^q+qq) leads to an even more complicated LCWF. 

To probe the parton content of the nucleon suitable models have to be invented to give explicit expressions for the 
LCWF amplitudes and exclusive processes have to be explored, where a large space-like momentum is transferred to 
an intact hadron. The framework for analyzing such processes was developed more than thirty years ago investigating 

■ elastic and inelastic form factors 0, H, Q and relies on perturbative QCD (pQCD). According to the factorization 
theorem, the scattering amplitudes can be expressed as convolutions of the (process-dependent but perturbatively 

■ calculable) hard kernel of the process and the nonperturbative (process independent) contribution describing the 
hadrons that take part in the reaction. In the case of form factors this contribution is represented by distribution 
amplitudes (DAs) @, H| that describe the hadron structure in parton configurations at small transverse separation. 
In the nucleon case, DAs probe that part of the nucleon state with orbital angular momentum L z — and at leading 
twist they involve only two of the six amplitudes entering the valence three-quark Fock component "0(3(3) ■ 

The properties of DAs were first studied using the method of QCD sum rules developed in Ref . § . This method gives 
the possibility to calculate the values of DA moments in terms of suitable sum rules. Therefore, knowing the behavior 
of the first few moments one can reconstruct the main properties of DAs as originally shown in Refs. [Tol. [ill. H3. fl3l [l4| . 
Although some work is available for other baryons (see, e.g., [l4l. HH [T 5 [l 7l. [Lsl. ITg| ) . the existing investigations were 
mainly focused on the nucleon DAs (see [ll| for an early review and [2 Oil for a more recent one) . Estimates of the 
nucleon DAs based on QCD sum rules can be found in Refs. [H, [H, [gj, [H, IH, [H, [H, HO, [3l|, H3L_ The nucleon 
DAs were systematically studied in Ref. [33| up to twist six and related to the nucleon form factors[34l. l35l , |3(| and 
the N — > A transition at intermediate values of the momentum transfer [3lj using light-cone sum rules. A variety of 
model calculations [37], [H, H^, H(| have also been derived from dynamical or phenomenological Ansatze for the nucleon 
wave function in order to describe the intermediate/low Q 2 region where the nonperturbative features of QCD are 
significant. Valuable additional information is also provided by lattice QCD [4ll. [42l. [43l. l44l. |45| . 

In other processes like deeply virtual Compton scattering (DVCS) or hard exclusive meson production the concept 
of generalized parton distributions (GPDs) [4(| has proven to be useful (for reviews, see J47[)- GPDs have been 
introduced as universal nonpertubative objects describing the hadron structure in terms of nondiagonal hadronic 
matrix elements of bilocal products of the light-front quark and gluon field operators. Their crossed version defines 
the generalized distribution amplitudes (GDAs) that describe the hadronization of a quark-antiquark or gluon pair 
in a pair of hadrons, e.g. a pair of 7r mesons, 7*7 — * 7T7t [48[. Other matrix elements can be defined generalizing the 
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concept of GPDs for non-diagonal transitions [49|, [50] and describing the transition amplitude between two hadrons 
or between a hadron and a real photon, thus called transition distribution amplitudes (TDAs). 

Recentl y, a ttention has been drawn to TDAs under the assumption that the factorization theorems for exclusive 
processes [5l| also apply to reaction mechanisms like proton-antiproton annihilation into two photons, pp — > 7*7, in 
the near forward region and large virtual photon invariant mass Q [52j or into a pion and a high-Q 2 lepton pair in 
the forward region, pp — > 7*7r — * l + 1~tt j53l |. exclusive meson pair production in 7*7 scattering at small momentum 
transfer 5J| , DVCS on a proton target in the backward region [55j , or hard exclusive electroproduction of a pion in 
the backward region, 7* A — » N'tt [56[. Within the factorization scheme, the hard and soft subprocesses decouple 
in the amplitude for these reactions, the soft part being a universal nonperturbative object describing the transition 
from a hadron to a real photon, or a proton to a pion. 

Depending on the values of the Mandelstam variables s and t in 7*7 scattering, a dual factorization mechanism 
has been identified in Ref. [57| describing the fusion of a real photon and a highly virtual and longitudinally polarized 
photon. One mechanism takes place when s <C Q 2 , while t is of the order of Q 2 , and involves twist-three GDAs, 
whereas the other one occurs for f < Q 2 and s ~ Q 2 and employs leading-twist 7 — » tt TDAs. Such TDAs have 
recently been studied in the Nambu- Jona-Lasinio model [{38| . The 7 — * tt~ TDAs are connected to the tt + — > 7 TDAs 
through CPT symmetry l59ll. The vector and axial- vector 71—5-7 TDAs have been analyzed in a quark model [6fj| , in 
the spectral quark model |6l| and in the Nambu- Jona-Lasinio model [62| . 

When studying the nucleon structure the N — ► tt TDAs are particularly interesting because they directly probe the 
three-quark plus sea qq pair component "0(3o±(?S-j n ©■ The possibility to extract experimental information on 
the N — > tt TDAs has been studied in Refs. [53.|63| for the pN — > 7*7r reaction in the kinematical regime accessible by 
GSI-FAIR [64[ and in Ref. [|36| for the 7* A — > N'tt reaction in the kinematical conditions of JLab. In such pioneering 
works the TDAs were predicted on the basis of the soft-pion theorems [65| which allow to calculate three out of the 
eight independent leading-twist TDAs in terms of the proton DAs. Predictions for the TDAs in the soft-pion limit 
were also obtained in Ref. [66| . However, it is desirable to extend these analyses to a more general framework, using 
as input different model calculations and also going beyond the kinematical soft-pion limit just because the A — > tt 
TDAs represent a new tool to access information on the Fock-space components of the nucleon wave function beyond 
the valence-quark contribution. Furthermore, in the impact parameter representation the A — > tt TDAs map out the 
transverse location of the small-size core and the meson cloud inside the proton 67 1. 

Being nonperturbative quantities, DAs, GPDs, GDAs and TDAs cannot be calculated from first principles, but have 
to be described by models or derived within lattice QCD. In this paper we are concerned with nucleon DAs and A — > tt 
TDAs within a phenomenological model for the LCWFs of the nucleon based on the light-cone constituent quark model 
(CQM) that has successfully been applied to the calculation of generalized parton distributions [68|, IfjjJ [7l[ , and 
transverse momentum dependent parton distributions [72l . [73j taking into account the full decomposition of the three- 
quark Fock component of the nucleon state. In order to derive expressions for the A — > tt TDAs we have to implement 
the model in order to include Fock components with a sea qq pair. This will be done assuming the nucleon to consist 
of a bare three-quark object surrounded by a meson cloud along the lines that were already followed in the calculation 
of the nucleon GPDs [74[ and electroweak form factors [75[. This will allow us to give a first estimate of the A — > tt 
TDAs for future applications. 

The paper is organized as follows. The nucleon DAs, whose definition and properties are recalled in Sec. [TTJ are 
explicitly derived in the light-cone CQM in Sec. IIIII and numerically computed in Sec. HVl The A — > tt TDAs are 
derived within the meson-cloud model in Sec. El and some results are presented in Sec. IVIl Concluding remarks 
are given in the final Section. The spin components required by the model for the baryon LCWFs are listed in the 
Appendix. 



II. NUCLEON DISTRIBUTION AMPLITUDES 



In this Section we recall some important definitions and properties of the nucleon DAs. 

In coordinate space, the proton DAs are derived from the following proton-to-vacuum matrix elements of trilocal 
operators built of quarks and gluon fields [l2], 0, [H[ 

(0|e yfe ?4 (zm)[zi; Zoli'iu'p {z 2 n)[z 2 ; Zo]fjd* (z 3 n)[z 3 ; z }k> k \p(p p ,X)), (2) 

where \p(p p , A)) denotes the proton state with momentum p p (p 2 = M 2 ) and helicity A; u, d are the field operators 
for up and down quarks, respectively; the Greek letters a, (3 and 7 stand for Dirac indices, while the Latin letters i, j 
and k refer to color; n is an arbitrary light-like vector (n 2 = 0) and z% are real numbers that specify quark separation, 
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with J2i Zi = 1. In Eq. J2]) the gauge factors [zi m , Zq] render the matrix element gauge-invariant and are defined as 



"Pexp 



ig(zi~z ) dtn^A* 1 (n[tzi + (1 - t)z 



(3) 



where V indicates the path-ordering prescription. In the following we will work in the light-cone gauge A + = where 
the gauge factors reduce to the identity. 

Taking into account Lorcntz covariance, spin, and parity conservation of the nucleon, the most general decomposition 
of the matrix element in Eq. ([2]) involves 24 invariant functions [33| . To the leading twist-three accuracy only three 
of them are relevant, and are given by the Lorentz invariant (scalar) functions of positive parity V(= vector), A(= 
axial- vector), and T(= tensor): 



u a( z i n ) uJ [}( z 2 n )dj(z3n)\p(pp, A)) = i/jv {fC) al3 {'~i h N + )~ f V(z t n ■ p) + {f lb C) al3 {N + ) 1 A{z i n ■ p) 

+ (<Jp„C) a f 3 (-f tl j 5 N + ) 1 T(z,7i -p) 



(4) 



where cr'" y = ^[7^,7^], cr PA1 is a shorthand notation for p v a v ^, C is the charge conjugation matrix and N + is the 
light-cone "good" or "large" component of the nucleon spinor N. In Eq. ((4]) the proton decay constant Jn is a 
dimensional quantity representing the value of the nucleon distribution amplitude at the origin of the configuration 
space [l(| EH, EH- Furthermore, we introduced a second light-like vector p^ such that 2p ■ n — 1. In particular, we 
make the following choice for the vectors p M and n M : 



p" = %l, 0,0,1), 



1 



(1,0,0,-1). 



It is convenient to define the functions V, A and T in momentum space 



V(x 1 ,x 2 ,x 3 ) = (n-p) 3 I Y[ ]Y)3 V ( Zin ' p > z z n 'P' Z3n-p)exp \i^2x k z k {n ■ p) 

■> 1 = 1 ^ n ' k=l 



(5) 



(G) 



and similarly for A and T. The variables x% conjugate to the light-cone positions of the quark operators in (U) are 
collincar momentum fractions of the proton longitudinal momentum carried by each quark in the infinite momentum 
frame, with < Xi < 1 and xi = 1 by momentum conservation. Accordingly, the scalar functions V(xi), A(xi) 

and T(xi) are DAs describing the longitudinal momentum distributions in the nucleon at a fixed scale p 2 that is not 
explicitly indicated if not necessary. 

Because of permutation symmetry between the two up quarks, the functions V and T are symmetric and A 
antisymmetric in their first two arguments. In addition, the requirement that the three quarks have to be coupled to 
give an isospin i state (the nucleon), yields the relations 



2T(1,2,3) = $(1,3, 2) + $(2, 3,1), 
$(1,2,3) = F(l,2,3)-A(l,2,3), 



(7) 

(8) 



which allow to express the proton DAs in terms of a single independent scalar function <f> with mixed symmetry. 

Introducing quark fields with definite chirality and denoting the Fourier transform of the matrix element on the 
left-hand side of Eq. ((4} by D ~, the three DAs can be obtained as 



V 
A 



1 



T = — 



(2)V4(p+)3/2 /Ar 

1 

- 2 l/ Hp + )3 /2 f N 
1 



(^12,1+ ^kl), 



^12,1 



4l,t 



D 



11.2' 



2V4(p+)3/2 f N 

where the |, J, arrows denote the up and down helicity of the proton, respectively. Thus, 



(9) 
(10) 

(11) 



(2)V4 (p +)3/2 /iv 



D 



12,1- 



(12) 
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Eq. (j4]) is equivalent to writing the three-quark uud component of the proton state with positive helicity in the 
infinite momentum frame as 1^, ID, EH 1 



J_fN 

V + A 



dx 



V -A 



u^Ch^i^dHh)} 



or in a more compact way as 

\p(P P ^)) uud 
where the integration measure is defined as 

da; 



J_f_N_ 

V3 4 





dx 




/■ 




10 




3 



«l(fci)«T(A 2 )dT(fc 3 )) _ T| U T(fc 1 ) u T(fc 2 ) rf l(fc 3 )) 



N 



> N 

n 



dx; 



N 



In Eqs. (JT3J) and flU 



ijk 



(13) 



(14) 



(15) 



(16) 



where b^ d (k, A) are the creation operators of free it and d quarks with momentum k = (k + , kj_) (fc + = (fc° + fc 3 )/\/2 
and being the plus and transverse momentum components, respectively), helicity A and color c (see also Eq. (|30[) 
below) . 

The corresponding neutron state is obtained from Q14p by interchanging u and d, with an overall change of sign. 

The matrix elements D^p , and ultimately the DAs, are directly linked to the L z = component of the valence- 
quark wave function of the nucleon, by integrating out the transverse momenta of the constituent quarks (l3l. [l4l. l37l. |39| 
(see Eqs. {32]) and 103} below). 

The DAs for the nucleon are well known in two limits [T2| . The first is the static SU(6) symmetric quark model, 
where the variables Xi take on only discrete values and the distribution amplitude is totally symmetric: 



$NR = S(xi - ±)5(x 2 - ^)S(x 3 - §). 
The second is the asymptotic regime of sufficiently high Q 2 where <J> takes the form 

$as = 120xiX 2 X3, 



(17) 
(18) 



which is totally symmetric under quark exchange and has the flavor-spin structure assumed in the SU(6)-symmetric 
quark model. In the asymptotic limit (Q 2 — > oo), A becomes negligible because of the Pauli principle, and V and T 
become totally symmetric under particle exchange, i.e. V,T — ► $as H3- 

Both limits (fl~7|) and (JTHJ) arc conflicting with experiment. In the first case one obtains wrong results for the neutron 
and proton magnetic form factors (G 7 ^ > 0, G P M < 0), in the second case > and G v M jG n M — > at large Q 2 fl2| . 

At intermediate values of Q 2 the nucleon DAs turn out to be quite different from their nonrelativistic and asymptotic 
limits. The analysis takes advantage of moments of the DA $ defined as 



^• m ' n) = f {dx}x[x™x^<t>( Xl ,x 2 ,x 3 ), 

4>N J 



(19) 



where [dx] = dxida;2dx3<5(l — x\ — x-2 — £3), and </>at is a normalization constant, which is chosen such that 0(°' o >°) = 1. 
Longitudinal momentum conservation {x\ + x% + X3 = 1) imposes the following condition 



(n 1 ,n 2 ,n 3 ) 



i(rii+l, n 2 ,n 3 ) 



j(m, n 2 +l,n 3 ) _|_ An-L,n 2 , n 3 +l) 



(20) 



1 Note that here the uud component is normalized as 

uud {p,\\p',\') uud = |2(27r)V%' + -P + )<5 (2) (P± -P±)<5a'A- 
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Thus, not all the moments at a given order M = ni + n<i + are linearly independent. 

The DA moments can be expressed in terms of matrix elements of suitable local operators entering appropriate 
sum rules 0, 0, llBl following the lines of the method of QCD sum rules developed in Ref. @. 

The nucleon DA obeys a renormalization-group equation which requires that ${xi,Q 2 ) is only logarithmically 
dependent on the momentum transfer scale Q 2 @. Following Refs. [13,123, the scale dependence of the nucleon DA 
can be cast in the form 



n=0 



log(Q 2 /A* CD ) 
log(^ 2 /A^ CD ) 



(21) 



where <& n (xi) are orthogonal eigenfunctions of the nucleon evolution equations, orthonormalized within a basis of 
Appell polynomials, and 7„ are the anomalous dimensions listed in Ref. In Eq. (|21[) the projection coefficients 
B n (fi 2 ) encode the non-perturbative input of the bound state dynamics at the factorization (renormalization) scale 
/i 2 . Using the explicit expression for the eigenfunctions $„ in terms of Appel polynomials they can be expressed as 
linear combinations of DA moments, i.e. 



120 



(22) 



where the coefficients a^- and the normalization constant N n have been calculated up to order M = i + j = 9 in 
Refs. P,[Z1- 

An alternative expansion of DAs is possible in terms of contributions of operators with a given conformal spin 

EE 

[l8l. |33[ |36| . This is convenient since operators with different spin do not mix under renormalization in one loop, and 
only operators with the same spin can be related by equations of motion so that the truncation of the conformal 
spin expansion at a certain order produces a self-consistent approximation. For example, at leading twist with the 
minimum possible conformal spin, <E> reduces to <&as , whereas its conformal expansion to the next-to-leading conformal 
spin accuracy reads [HI, Hfjj] 



${ Xi , M 2 ) - $as(^)<^ 2 ) 1 + <^(M 2 )(zi - x 2 ) + ^(M 2 )(l - 3x 3 ) 



(23) 



where 



IN- 



21 r 



0(1,0,0) _ 0(0,1,0) 



4.T 



1,0,1) 



(24) 



Numerical estimates of the coefficients at fi 2 = 1 GeV 2 available from QCD sum rules [H, EH, HH, H3| give 



f N = (5.3 ± 0.5) x 10" 3 GeV 2 



= 4.0 ± 1.5, 



1+ = 1.1 ±0.3. 



(25) 



An approximately 40% lower value of /jv has been determined recently in lattice calculations [ij] : fjq = 3.234(63) (86) x 
10~ 3 GeV 2 at [i 2 = 1 GeV 2 . In our numerical estimates the value in Eq. (f2"5|) will be used. 

The resulting DA exhibits a broad and rich structure that is reflected in an asymmetric distribution of the proton 
momentum between the three valence quarks (in the limit of infinite momentum). According to the QCD sum rule 
approach about 60% of the proton longitudinal momentum is carried by the up quark with its helicity in the same 
direction as that of the proton. The remaining up and down quarks, with combined helicity zero, are confined into the 
small- x region, each carrying about 20% of the total longitudinal momentum. This asymmetry is a common feature 
of all octet baryons [T3 |. A somewhat smaller asymmetry of the helicity amplitude tit for the uud configuration was 
found in Ref. [23j re-evaluating the momentum sum rules. This suggests the possibility of spin-zero diquark clustering 
in the nucleon wave function as a manifestation of the attractive QCD force produced by gluon exchange that is just 
strongest in the spin-zero quark-quark state [39l |. However, lattice QCD calculations of the first two moments were 
unable to confirm this asymmetric behavior of the nucleon DAs ■ 

The problem with QCD sum rules is that the moment sum rules are not stringent enough to fix the shape of 
the nucleon DAs uniquely 0, [12| ■ With increasing order of expansion the oscillations become stronger, and small 
variations of the moments may lead to a completely artificial behavior (25[. Such an extreme sensitivity of the 
expansion coefficients indicates that the moments do not give a convergent expansion. Actually, there is an infinite 
number of possible solutions which satisfy the moment sum rules, but differ dramatically in their shape. Thus, the 
predicted observables, like form factors, may result quite different. To determine the possible variation of DAs allowed 
by moment sum rules and to reconcile the constraints from moment sum rules with data, the heterotic solution was 
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proposed in Refs. [27], [28[ combining Q 2 evolution equations in pQCD, QCD sum rules and phenomenology. Allowing 
some flexibility to the expansion coefficients (|22[) . while keeping them as close as possible to the sum-rule requirements, 
a good agreement with high-Q 2 data on the proton magnetic form factor was achieved. 

Alternative phenomenological approaches take advantage of constraints on the three-quark component of the nucleon 
wave function imposed by data. For example, requiring that with the same nucleon wave function one reproduces 
the proton form factors, the phenomenological valence quark distribution as well as the Jip — ► pp decay width, the 
authors of Ref. [40| assumed a wave function with a small hard factorizing part depending on the longitudinal momenta 
and described by a nucleon DA and a large soft nonfactorizing contribution depending on Xi and kjj_ solely in the 
combination k 2 j_/xi with a Gaussian fall-off with ka_. All requirements are met with a model wave function depending 
on only two parameters, namely the proton decay constant /jy and the size parameter of the transverse momentum 
dependence. This gives a DA that is much less asymmetric than that derived from QCD sum rules, rather resembling 
the asymptotic DA, but with the position of the only maximum somewhat shifted. Actually, the asymmetry of the 
leading-twist amplitude constraint by phenomenology [HI, [36|, H(| or calculated on the lattice [HI, 0, EH is much 
smaller than in QCD sum rule calculations (see Table U discussed in Sec. IIV[) . 



III. NUCLEON DISTRIBUTION AMPLITUDES AND LIGHT-CONE WAVE FUNCTIONS 



In this Section we derive explicit expression for the matrix elements D *n in terms of LCWFs. To this aim, we 
first introduce the Fourier expansion in momentum space of the quark field operator of flavour q and colour c [l|, Q 

q (zn ,zn ±) = y__e(*+) 

x *^^{bg(k, X)u+(k, A) exp(— ik + zn~ + ik± ■ zn±) 
x 

+ d' q c (k, X)v+(k, A) exp(+jfc + zn~ — ikj_ ■ zni)}, (26) 

where the b and d' operators respectively annihilate the "good" component of the quarks fields and create the "good" 
component of the antifields fulfilling the anticommutations relations 



{^',(fc',A'),^ c (fc,A)} = {^',(fc',A'),4 c (fc,A)} 

= 16^ 3 fc+<5(fc'+ - fc+),5( 2 )(k' ± - k x )V<A'A<W 

In Eq. (|26p . u+(fc, A) and v + (k, A) are the light-cone spinors of the quark and antiquark, respectively. 
The three-quark Fock component of the light-front proton state is given by 



(27) 



\p(p P ,X))= E / ^1 [d'kJa^'^ae.k^jA,,^}^,^)!]!^;^,?^), 



(28) 



where \&?'^({£j,ktj_; Xi, Ti}i— 1,2,3) is the momentum LCWF which gives the probability amplitude for finding in the 
nucleon three quarks with momenta (CiPp , Pj_l — k^j_ + CiPp)i an d spin and isospin variables A^ and Ti, respectively. 
The proton state is normalized as 



(p,X\ P ',X) = 2{2nfp+5{p l+ -p+)SW(p' ± -p ± ), 
and the three-quark state is defined as 



c ijk 



ni9 Al ;^) = ^^(fci,Ai)6^(fc 2) A 2 )6t fe (fc3,A 3 )|0). 



i=l 



(29) 



(30) 



In Eq. (|28p and in the following formulas, the integration measures are defined by (|15p and 



N 



d 2 k,; 



iV 



(31) 
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Inserting the momentum-space expansion ((26)) of the quark fields and the proton Fock-statc (|28[) in Eq. and 
using the anticommutation relations for the quark creation and annihilation operators, one obtains for the matrix 



elements 



a/3,7 



-24- 



1 



-_u +a (xip+, \i)u+0{x2Pp , A 2 )u +7 (x 3 p+, A 3 ) 



pAf] I {si.ku.; Ai, 1/2}, {x 2 , k 2± ; A 2 , 1/2}, {x 3 , k 3± ; A 3 , -1/2} 



(32) 



The light-cone spinors in Eq. 



are explicitly given by 



1 ZiPp 

V2 



and 



-(^tPp,J.) 



1 £iPp_ 
V2 



V 



'J 



1,2,3. 



(33) 



As a consequence, the Dirac component of the quark spinor selects the quark-spin configuration, with a = 1,3 
(a = 2, 4) corresponding to helicity A =f (A =|) of the quark. In particular, from Eqs. (f9l)- (fTTj) we see that for a proton 
with helicity | the DAs are obtained from the wave- function component with total quark helicity ^, corresponding 
to the projection onto the partial wave with orbital angular momentum L z = 0, and to the three spin configurations 
(TTI)j (TIT)? an d (ITT)- Furthermore, the quantum numbers for the quark isospin in the LCWF in (|32[) correspond 
to the isospin projection in the uud configuration. 



IV. NUCLEON DISTRIBUTION AMPLITUDES IN A LIGHT-CONE QUARK MODEL 

In this Section we estimate the DAs using a phenomenological model [H, [H, [73, [ill, [71, [73| with LCWFs built 
in such a way as to satisfy Poincare covariance and to be eigenstates of the total angular momentum operator in 
the light-front dynamics. These properties can be fulfilled by constructing the wave function as the product of a 
momentum-dependent wave function tp({xi, k^j_}) which is spherically symmetric and invariant under permutations, 
and a spin and isospin wave function which is uniquely determined by symmetry requirements and invariant under 
permutations, i.e. 

^({xuk^Xun}) = ^({x,,k 4 x}) 

>< E <i;(^/(fci))^S;(^/(fc 2 ))< 3 2 ; 3 (i? c /(fc3)) 

A«lA'2A'3 

X $*l(Ml)M2,M3,Tl, T 2 ,r 3 ), (34) 

where D 1 ^ (R c f(k)) arc matrix elements of the Melosh rotation R c f [7^|, which converts the rest-frame spin of the 
quarks into light-cone spins. They are explicitly given by 



Dx*(Rcf(k) = (\\R c f(k)\ri 

where m is the quark mass and M is the mass of the non-interacting three quark system. In Eq. (|34[) the spin and 
isospin quantum number of the quarks are coupled by the SU(6) symmetric function defined as 

*Xr(Mi) M2,M3,Ti, r 2 ,r 3 ) = -^=[$i (^1,^2,^3)^1 (n,T2,T 3 ) + $i (mi,M2,M3)*i (n , T 2 , r 3 ) 1 , (36) 

V2 L 2 A 2 r 2 A 2 T J 

where 

*j"= E^/ 2 '^ ;1 /2,M2ki2,M Jl2 )(J 12 ,M Jl2 ;l/2, M3 |J,A). (37) 

M Jl2 

In the case of the proton DAs, we need the uud isospin projection of the nucleon wave function, which corresponds to 
^({^k^AJ,}^}) = -4^({a; l ,k 8± })SP(Ai,A 2 ,A 3 )l> 1 11 (ii -±), (38) 
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with the isospin coefficient 1(5151 — ■§■) — A/'f an< ^ the s P m dependent part given by 

22 

s5(Ai,a 2 ,a 3 ) = Yl <A;(^/(fci))<i;(i? c /(fc2))< 2 A ;(i? c /(fc 3 ))$ i i (mi,m2, 



(39) 



Table I: Results for the moments 0"' m ' n ' with I + m + n < 3 of the proton DA in different model calculations: COZ from 
Ref. H; KS from Ref. [H|, SB from Ref. DF from Ref. HH, BK from Ref. [U, and PPB from the present model. 



{1, 


m 


n) 


COZ 


KS 


SB 


DF 


BK 


PPB 


(0 





0) 


1 




1 


1 


1 


1 


1 


(1 





0) 


0.54 - 


0.62 


0.46 - 0.59 


0.572 


0.582 


0.381 


0.346 


(0 


1 


0) 


0.18 - 


0.20 


0.18 - 0.21 


0.184 


0.213 


0.309 


0.331 


(0 





1) 


0.20 - 


0.25 


0.22 - 0.26 


0.244 


0.207 


0.309 


0.323 


(2 





0) 


0.32 - 


0.42 


0.27-0.37 


0.338 


0.367 


0.179 


0.151 


(0 


2 


0) 


0.065 - 


0.088 


0.08 - 0.09 


0.066 


0.085 


0.125 


0.141 


(0 





2) 


0.09 - 


0.12 


0.10 - 0.12 


0.170 


0.083 


0.125 


0.136 


(1 


1 


0) 


0.08 - 


0.10 


0.08 - 0.10 


0.139 


0.108 


0.101 


0.099 


(1 





1) 


0.09 - 


0.11 


0.09 - 0.11 


0.096 


0.106 


0.101 


0.096 


(0 


1 


1) 


-0.03 - 


- 0.03 


unreliable 


-0.021 


0.018 


0.083 


0.091 


(3 





0) 


0.21 - 


0.25 




0.21 


0.249 


0.095 


0.076 


(0 


3 


0) 


0.028 - 


- 0.04 




0.039 


0.041 


0.059 


0.070 


(0 





3) 


0.048 - 


0.056 




0.139 


0.040 


0.059 


0.067 


(2 


1 


0) 


0.041 - 


0.049 




0.079 


0.060 


0.042 


0.038 


(2 





1) 


0.044 - 


0.055 




0.049 


0.059 


0.042 


0.037 


(1 


2 


0) 


0.027 - 


0.037 




0.050 


0.040 


0.036 


0.037 


(1 





2) 


0.037 - 


0.0434 




0.037 


0.039 


0.036 


0.035 


(0 


2 


1) 


-0.004 - 


- 0.007 




-0.023 


0.004 


0.030 


0.034 


(0 


1 


2) 


-0.005 - 


- 0.008 




-0.007 


0.005 


0.030 


0.033 



. [d 2 kJ 3 V({^,k I± }) s?(U,T) + 3?a,T,T) 

In J L 1 ■ 

J [d 2 k ± ] 3 ^({^, k^}) [sf a, t, T) - s?(t4. T) 

4^3 



Inserting the wave function p8|) in Eq. (|32|) . we hnd for the nucleon DAs the final results 

v A ^ 

A 
T 
4> 



TN 



8\/3 



IN 



d 2 kJ 3 ^({^,k i± })Hf(T,T,|), 
[d 2 kJ 3 ^({x,,k^})^(T,U), 



(40) 

(41) 
(42) 
(43) 



where the explicit expressions for the function arc given in the Appendix. These results confirm that by integrating 
out the transverse- momentum dependence of the nucleon wave function, DAs are determined by its L z = component. 
In particular, in the present SU(6) symmetric model they involve only one of the two independent light-cone amplitudes 
parametrizing the S-wave component of the LCWF. In order to probe also the other light-cone amplitude one should 
consider a more general framework with mixed-symmetry terms [80| . 

In the following the results of Eqs. (|40|) -(|43ll arc applied to a specific CQM taking the form of the momentum wave 
function from Ref. [8L|, i.e. 



i>({ Xi ,k i± }) = 2(2tt) 2 



1 LU1LO2CO3 
M X1X2X3 



1/2 



N' 



(M 2 +/32)T' 



(44) 
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where u>i is the free-quark energy and N' is a normalization factor such that J d[x]3|^({xi}, {&i_i_})| 2 = 1. In Eq. I|44[). 
the scale /?, the parameter 7 for the power-law behaviour, and the quark mass m are taken from Ref. (8l| , i.e. (3 = 0.607 
GeV, 7 = 3.4 and m = 0.267 GcV. According to the analysis of Ref. [12] these values lead to a very good description 
of many baryonic properties. 

In Fig. [T] the model results for the proton distribution amplitude $ are shown. The resulting shape is quite similar 
to that of the symmetric asymptotic DA in Eq. (|18|) . 








Figure 1: The proton distribution amplitude. 

In Table 1 we list the 0('> m >") moments up to the order M = I + m + n < 3 m comparison with different model 
calculations. At this order there are 20 moments out of which only 10 are independent. Despite the fact that $ is 
normalized, one has to keep in mind that it is a distribution amplitude, not a probability density. Thus its moments 
cannot strictly be interpreted as mean values [27| . However, the first moments of 4> provide an indication on how 
the longitudinal momentum of the proton is partitioned among the valence quarks. Since tp^ 1 ' ^ sa 0.6 in the model 
of Ref. [34||, this is the reason why it is claimed that with QCD sum rules about 60% of the proton longitudinal 
momentum is carried by the up quark with helicity parallel to that of the proton. From Table [J we see that a quite 
different result is obtained in the light-cone quark model, a roughly symmetric result similar to the one considered in 
Ref. j43| . In fact, while starting from quite different assumptions, the LCWF used here and the model wave function 
of Ref. [4(| have some common features. Besides the longitudinal momentum dependence of L z = part of the 
present LCWF giving rise to an almost symmetric DA, our model resembles the Gaussian fall-off as function of the 
quark transverse momentum at large x of the wave function of Ref. and both models are able, for example, to 
reproduce the nucleon form factors. As a matter of fact, the transverse-momentum dependence of the present-model 
LCWF has been studied in Ref. [73[ , showing that the ratio of the squared mean transverse momentum and the mean 
square transverse momentum of transverse momentum dependent parton distributions agrees within 10% accuracy 
with the results obtained assuming a factorized gaussian form for the nucleon wave function. 

Using these results, we can study the behaviour of proton DA calculated in our model after evolution from the 
initial low-scale of the model to higher scales. The initial scale, corresponding to the results shown in Table [IJ has been 
fixed evolving back unpolarized data, until the valence distributions matches the condition that the second moment, 
i.e. the momentum fraction carried by the valence quarks, is equal to one. Using LO evolution equations, we find 
Qq = 0.079 GeV 2 , with Aqcd = 0.232 GeV [6{|. Although there is no rigorous relation between the QCD quarks 
and the constituent quarks, and a more fundamental description of the transition from soft to hard reg imes would be 
very helpful, this strategy reflects the present state of the art for quark model calculations [H, [H, |85| |. and has been 
validated with a fair comparison to experiments (see, for example, Ref. [z3). The results for the first moments <^> m,ri ) 
with I + m + n < 2 after evolution to Q 2 = 1 GeV 2 are shown in Table [Til Comparison with the analogous results 
from lattice estimates [45[ is quite nice. Evolution has only a very small effect: a further evolution to Q 2 = 4 GeV 2 
would only hardly modify the last digit of our result. This is in agreement with the logarithmic scale dependence 
predicted by Eq. (j2"Tj) and the fact that already at the input scale the behaviour of our DA approaches that of the 
asymptotic DA. 
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Table II: The moments 0('' m,n ) of the proton DA at different scales. The asymptotic values (AS) of Eq. (| 1 8fl (second column) 
are compared with the lattice results (LAT) [3lJ at the scale Q 2 — 1 GeV 2 (third column) and with the corresponding results 
from the present model calculation (PPB), after evolution from the initial scale Q 2 , — 0.079 GeV 2 to Q 2 = 1 GeV 2 (last 
column) . 



(I, m, n) 



AS 



LAT 



PPB 



(o 





o) 




1 


1 


1 


(1 





o) 


1 

3 


~ 0.333 


0.3999(37)(139) 


0.340 


(0 


1 


o) 


1 
3 


~ 0.333 


0.2986(11)(52) 


0.335 


(0 





1) 


1 

:) 


~ 0.333 


0.3015(32) (106) 


0.326 


(2 





0) 


1 
7 


~ 0.143 


0.1816(64) (212) 


0.147 


(0 


2 


0) 


1 

7 


~ 0.143 


0.1281(32) (106) 


0.144 


(0 





2) 


1 
7 


~ 0.143 


0.1311(113)(382) 


0.137 


(1 


1 


0) 


2 

121 


~ 0.095 


0.1092(67) (219) 


0.098 


(1 





1) 


2 

21 


~ 0.095 


0.1091(41) (152) 


0.095 


(0 


1 


1) 


2 

21 


~ 0.095 


0.0613(89)(319) 


0.093 



TRANSITION DISTRIBUTION AMPLITUDES IN A MESON-CLOUD MODEL 



The general matrix clement describing the transition from a nucleoli to a meson state reads |5£ 
(ir\e l: > k q l a (zin)[zi; zo^'iq 3 jp{z 2 n)[z2; z ]j> jq* (z 3 n)[z 3 ; z ] k > k\N) , 



(45) 



where the Wilson lines [zf, zq] are defined as in Eq. ([3]). In the following they will be neglected by assuming to work 
in the light-cone gauge A + = 0. The spinorial and Lorentz decomposition of the matrix element (|45[) follows the same 
line as in the case of the baryon DAs. In particular, for the p — ► ir° transition the leading-twist TDAs can be defined 
as 



4^((7rV)|e <J '*«a(*i»)«M*2n)^(« 3 n)b(ft, I A)) 



Jn 



Vr {fC) aP {N + ) 1 + (^ 5 C)^( 7 5 iV+) 7 

+ Tf\a w C) a p{ 1 ^N + ) 1 + M-% pw ° (fC) ct0 (4L x N+) 1 
+ M- l A% °(^ 5 C) a0 (r^ ± N+) 1 + M~ 1 Tr°(<J pA± CU(N + ) 1 
+ M- 1 ^ {o m C) aP {o^^N+) 1 + M- 2 Tr° (a pA± C) a p(4L ± N + )~ 



(46) 



where the symbol T represents the Fourier transform (like Eq. ([5])) and f n is the pion decay constant (/ w = \f ! !F. K = 131 
MeV) . In a reference frame with the z-axis along the direction of the proton momentum, the pion momentum p n and 
the proton momentum p p have the following Sudakov decomposition 



M 2 

p P = (1 + 0+ Y+~£ n ' 

P, = (l-0p+ ml + Al n 



where A is the four-momentum transfer, 

A = p n -p P = -2£p 



i-e 



A 2 



Ax, 



i-e 



M 2 

1+c 



(47) 
(48) 

(49) 
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and £ is the skewness variable describing the loss of plus momentum of the initial hadron in the proton-to-meson 
transition, i.e. 



A • n A+ . , „ 1 , 

C = = n with P = —\Vn + P-k)- 

^ 2P-n 2P+ 2 W w 



(50) 



The TDAs are dimensionless functions and depend on (xi, £, A 2 ), where the fraction of quark plus momentum Xi have 
support in [— 1 + £, 1 + £] and 



A 



-2C 



M 2 



Ai. 



(51) 



Restricting ourselves to the case £ > 0, momentum conservation requires Xj = 2£. The fields with positive 
momentum fractions. x, > 0, describe creation of quarks, whereas those with negative fractions, a;, < 0, the absorption 
of antiquarks. This leads to define three distinct kincmatical regions: the ERBL region for x, > 0, and two DGLAP 
regions when x\ > 0, X2 > 0, X3 < 0, or xi > 0, X2 < 0, X3 < 0. The names derive from the evolution equations which 
controls the scale dependence of the TDAs in the different regions. 

Denoting the matrix element in left-hand side of Eq. (|4"6")l by „, we can derive the eight TDAs in terms of the 
following linear combinations of matrix elements 



rpPTTO 



rppiTo 



1 



'2 1 /4yr+^(p+)3/2 /Ar 
1 u 



T 



■T 1 

^21,1 



'2 1 /4yi + ^(p+)3/2 /jv 

1 u 

'2i/4yrT^(p+ )3/ 2 /jv 

MAT 1 



I 1 12,1 ^21,1 J ' 



J 11.2 



Ai 



2 J 22,2 



/tt 



2 1 /4yTTe(P+) 3 / 2 



A/A" 



1 



ai 2i/4yr+c(^+) 3/2 /. 



iV 



1 



A^ 2 1 /4yr + ^(P+)3/2 f N 

M 1 



"T ^21,2 



^12,2 ^21,2 )) 



A + T^T _ A rp\ 

11,1 ^±^22.1 



-Tt 



'Ai 2 1 /4^/T+£(p+)3/2 f N 

,2M 2 (Ai) 2 1 



1 11,1 ' ^J. J 22,l 



./V 



(Ai) 2 2i/4^T+^(p+)3/2 ; 



N 



22,2: 



(52) 
(53) 
(54) 
(55) 
(56) 
(57) 
(58) 
(59) 



where A^ = A x ± iA y . 

In the following we focus on the study of the TDAs in the ERBL region, corresponding to probe the ijj3q+qq Fock- 
componcnt of the nuclcon wave function. The five-parton component of the nucleon state can be modeled using the 
meson-cloud model developed in Refs. [7J, [75[. The basic assumption of the model is that the physical nucleon is 
made of a bare nucleon dressed by the surrounding meson cloud so that the nucleon state is decomposed according 
to the meson-baryon Fock-state expansion as a superposition of a bare nucleon, formed by three valence quarks, and 
states containing virtual mesons with recoiling baryons. These baryon-meson subsystems are assumed to include 
configurations with the baryon being a nucleon or a A and the accompanying meson being a pion as well as a vector 
meson such as the p or the to. Being interested to the p — > 7r° TDAs, here we will consider the meson-baryon 
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components with a pion and a nucleon or a A, given by the following representation in the light-cone dynamics 

(» - E e.) f 3 ' (ki - 1 k '^) *(i - E f.) *' 2 ' fix) 

X E E 

A' Xi,Ti,Ci 

5 

X (1 - y, -k ± ; k^; A i; 7-;}; =4 , 5 ) l&Pp > + &Pp-L> A " T *' ^ ?>> ( 60 ) 

i=X 

where the LCWF of the baryon, S&f,, and the pion, $ 7r >[/] ) incorporate the Jacobian of the transformation from 
the intrinsic variables with respect to the hadron rest-frame ({Q, to the intrinsic variables with respect to the 

nucleon rest frame k'-j_}), i.e. 

§f; 1/1 ({&, k^ ± ; Ai, n} i= i,2,3) = \A7(6,k,k)*f/ [/1 ({0, A*, n} i= i,2, 3 ) 

= -^^'^({Ci.KtL;^, r < } i =l I 2,3) ) (61) 



^({^k^A;, -7-^=4,5) 



[/1 ({Ci, «U.; A*, Til^s) 

TT^r*"' 1/1 ({0, K iX ; A,, 7i} i= 4 l6 ). 

l - y 



The relations between the variables ({cci,kix}) and {£.;,k' i± }, arc given by: 
For i = 1, 2, 3, corresponding to the indices of the three quarks in the baryon, 



= K± - Ci k J 



For i = 4, 5, corresponding to the indices of the quark and antiquark in the pion, respectively: 

Ci 







(i-y)' 



= k- ± + Ci kj_ 



(62) 



(63) 



(64) 



In Eq. (|60p the function tf)^^' B7r \y,]i±) is the probability amplitude to find a physical nucleon with helicity A in 
a state consisting of a virtual baryon B = N, A and a virtual pion, with the baryon having helicity A', longitudinal 
momentum fraction y and transverse momentum k^ , and the pion having longitudinal momentum fraction 1 — y and 
transverse momentum — This splitting function can be calculated using time-ordered perturbation theory in the 
infinite momentum frame as explained in Ref. [86j , and have also been rederived and tabulated in Ref. [74| . 

For the pion state in the matrix element of Eq. (|46|) we consider the valence qq component given by 



Ti ,Ai ,Ci 



^)}= E / % [d 2 kx]2^ [/1 (U„k l x;A l ,r l } l=1 , 2 )n|g^;^+, Pl 



(65) 



where pj 



£,iPn and the qq state is defined as 



i=l 



^&ti(fc 1)Al)d « ( fc 2)A2) | }. 



(66) 



Using the expressions for the proton and pion state given in Eqs. (|60p and (|65[) . and the momentum-space expan- 
sion (|26|) of the quark fields, combined with the anticommutation relations for the quark creation and annihilation 
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operators, the final expression for the matrix elements ^ is given by 

/ i \ 3/2 j 

T a,/3 7 = - 24 ^7 ) 7 M+a(fci + , Ai)u +/3 (fc+, A 2 )w +7 (fc+, A3) 

x I [d 2 K ± } 3 *v 1/1 ({f|, kix; Ai, 1/2}{|, k 2± ; A 2 , 1/2}{|, k 3 x; A3, -1/2}) . 

(67) 

The light-cone spinors of the quarks are defined as in Eq. (|33[) . and depend on the longitudinal momenta = XiP + . 
The Dirac indices a, (3, and 7 fix the total quark hclicity of the baryon wave function, as explained in sect. [Til while 
the isospin quantum numbers in the baryon wave function correspond to the uud configuration. Eq. (|67p has a clear 
physical interpretation and allows to relate the nucleon-to-pion TDAs to the baryon distribution amplitudes in the 
(Bit) component of the nucleon weighted by the probability amplitude that the nucleon fluctuates in the corresponding 
(Bit) subsytem with the pion momentum matching the pion momentum in the final state. The momentum fraction 
of the quarks in the baryon LCWF are defined with respect to the longitudinal momentum of the baryon, i.e. 
K t /(VPp) = x i/{2£,)i while the integration over the transverse quark momenta corresponds to the projection of the 
baryon LCWF onto the zero orbital angular momentum component. In Eq. (|67p . the sum over the baryon states is 
restricted to the nucleon and the A, while the sum over the helicity A' of the baryon permits baryon-pion fluctuations 
which do not conserve the hclicity of the parent nucleon. 

In the case of the nucleon contribution, we model the proton LCWF as in Eq. (f3"8")) . with parameters 7 = 3.21, 
P = 0.489 GeV and m = 0.264 GeV from the fit of the valence and meson-cloud contribution to the electroweak 
nucleon form factors [75j]. The A is described as a state of isospin r = | obtained as a pure splin-flip excitation of 
the nucleon, with the same momentum-dependent wave function of the nucleon, i.e. 

v^- [/1 ({xJ,{fcJ,{AJ,{ MM 4) = ^({^,k J± })S^(A 1 ,A 2 ,A 3 )$ 1 31 (i,i,-i), (68) 

2 2 

with the isospin coefficient ^31(^,5,-^) = \ \ and the spin dependent part given by 

22 

E£(Ax,A 2) A 3 )= ^A*(^/(fci))^^ (69) 

The explicit expression of the functions for all the possible spin configurations of the three quarks in the A state 
is given in the Appendix. 

Finally, the splitting function 4>^' B ^ in Eq. (|6~7| is calculated as in Ref. [75^ 

Inserting the matrix elements p from Eq. (|6~?| into Eqs. (|52")) - (p)9")) we finally obtain the expressions of the TDAs 
in the meson-cloud model: 
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x {ct>][ N ^\y, -p w± ) [sf (1,1, T) + Hf (I, T, t)] + 0j^' A7r) (z/, -p„ x ) [sf (T, 4, T) + sf a, T, T)] } , (70) 



AT = -i^^-JpA /^^({^l) 



{t}i N ' N *\y, -p^) [sf (U, T) - Sf a, T, T)] + ^r ,A7r) (z/, -p w x) [sf (t, I, T) - sf (i, T, t)] } , (71) 
d 2 k ± ] 3 ^({^, k 4± }) {<f>\i N ' N%) (y, - Pw± )sf (t, T, I) 



■ 4V3 / (1 - Q A 
1 2£ V(l + 3 /iv 

+0;^ A7r) ( y ,- P7r± )sf au) + ^4^^ 



^ - ^v^S^ry [dk ^ ({a;i ' k - }) 

x {^^(y, -p, J [sf (t,|, I) + sf (|,T, I)] + 4r ,A7r) (y, -p„_l) [sf (T, 4, 1) + sf a, T, 4)] } , (73) 



Ar = ^/t^^T 1 ^ /[d 2 k ± ] 3 ^,k l± }) 



.4^3 / (1 - Q U MA" 
2£ V (1 + £) 3 /jv Ai 

x {^-^(y, -p, J [sf (t, I, I) sf (4, t, I)] + ^' A %, -p.±) [Sf (T, I, I) - sf (4, T, 4)] } , (74) 

= ^^^twj^ I [ d2k J3^, kii}) {-a: ^^(y, - Pir x)sf a, i,t) 
+a+ A %, - p . ± )s a (t, t, t) - a : ^^ A7r) (y, - P ^)sf aa, t)} , (75) 

?T = -^^JjTTW^ J [d 2 kJ ; #({x,,k 4± }){A I ^^( yi - P ^)sf (14 >T ) 

+A+ ^'^(y, - P .JS A (T, T, T) + Aj ^' A7r) (y, -ivJsf (U, T)} , (76) 



.4^3 /(1-0 / w 2M 2 (A] 



\2 



IT" - -^^^T#7^ "'(A1)^ 7 [d 2 kx]3^({^ ; k^})4 ( o JV!A7r) (y,-p 7r± )S A (|,|,l), (77) 

where the longitudinal momentum fraction in the argument of the splitting function is y = 2£/(l + £) and the spin | 
(— |) state of the A is indicated as -ft- (JJ.). 



VI. SOME RESULTS FOR THE TDAS 



As an example results are reported in Figs. [2] and [3] under kinematic conditions relevant in the case of hard exclusive 
electroproduction of a pion in the backward region [56| or in associated production of a pion and a high-Q 2 dilepton 
pair in pp annihilation [63] . Preliminary results for other kinematics have been presented in Ref. (88j . 

The vector and axial- vector TDAs in Fig. [2] exhibit the expected symmetric and antisymmetric behaviour under 
permutation of the two up quarks, respectively. The A contribution has the same shape as the proton contribution, 
with the same sign for A v ^ and , and opposite sign for Vf* and A!% . The relative contribution of the nucleon 
with respect to the A is always smaller in absolute value for the vector TDAs, being suppressed by a factor of about 
10 in the case of V/" r and by a factor of about 1.5 in the case of V 2 p7r . Viceversa, for the axial- vector TDAs one finds 
that the nucleon contribution to AP^ is smaller than the A contribution by a factor 3, while for ASf the weight of 
the nucleon contribution is three times larger than in the case of the A. This can be traced back both to the different 
spin structure and to the different splitting functions in Eqs. (|70|) - ([7T|) . and (|73|) - (|74p . In particular, V[ v and 
involve splitting functions without flip of the helicity of the parent nucleon, while VJ^* and are proportional to 
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splitting functions with helicity flip. In the explored kinematics, the (AT, Ntt) vertex without hclicity flip is suppressed 
by a factor 5 with respect to the (TV, An) interaction, while for the opposite case with hclicity flip the probability 
amplitude to have a (Nir) fluctuation in the nucleon is almost twice larger than for the (A71") subsystem. 

The tensor TDAs in Fig. [3] are symmetric under permutation of the two up quarks and, at variance with the other 
TDAs, involve spin configurations also with parallel helicities of all three quarks. Such spin configurations receive 
contribution only from the A with helicity J z = ±3/2 because of the projection of the baryon wave function in 
Eqs. (|72[) and l|75|) - (|77|) onto the zero orbital angular momentum component. In particular, one finds that in the 
explored kinematics the splitting function of the nucleon into a A with helicity J z =3/2 is 10 times bigger than in 
the case with hclicity J z = —3/2, and with opposite sign. Furthermore, these terms are multiplied by kinematical 
coefficients which modulate their relative contribution to the tensor TDAs in a quite different way. For example, 
the contribution of the A state with helicity J z =3/2 is quite small in the case of Xf', while it is enhanced by an 
additional factor of 2M 2 /Aj_ ~ 50 in T , J" r . This is the only non-vanishing contribution to T|" r . In Tf™ there is also 
the contribution from the helicity states J z = 1/2, with the same sign for the nucleon and the A, but larger by about 
a factor 10 in the case of A. In the case of Xf™ and the contribution from the A is large, and it is mainly given by 
the configuration with hclicity J z = 3/2. However, this contribution is reduced by that of the proton with a similar 
shape in the case of Tf 71 ^ whereas 2/3 of Tj™ are due to the proton and 1/3 to the A. 




Figure 2: The p — » n° transition distribution amplitudes Vf* (up left), (up right), V 2 p7r (down left), (down right) as 
function of (xi, X2, 2£ — x\ — x%) at fixed £ = 0.9 and A 2 = —0.1 GeV 2 . 



VII. CONCLUDING REMARKS 

In the light-cone description the nucleon state is decomposed in terms of A^-parton Fock states with coefficients 
representing the momentum light-cone wave function of the N partons. Since the constituent quark models work so 
well phenomenologically, in applications it is usually assumed that only the Fock components with a few partons have 
to be taken into account. One of such models has been studied in a series of papers to show that the parametrization of 
the LCWF up to five-parton components is already sufficient to account for the electroweak form factors [75[ and spin 
densities [7pJ of the nucleon, as well as the observed asymmetries due to transverse momentum dependence of parton 
distributions (zl, [73| , and to give a useful insight into the quark generalized parton distributions [68|, [69|, [7l|, [74[ . 
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As a further test of the model in this paper the nucleon distribution amplitudes and the nucleon-to-meson transition 
distribution amplitudes have been considered. At leading twist the nucleon DAs probe the three-quark content of 
the nucleoli state 1 with orbital angular momentum L , = and the N — tt TDAs probe the qq sea pair contribution 
responsible for the meson cloud surrounding the bare three-quark nucleon. 

Assuming SU(6) symmetry the shape of the calculated nucleon DA is similar to the asymptotic DA, with a roughly 
symmetric contribution of the three quarks. This contrasts the results from QCD sum rules that push towards highly 
asymmetric quark contributions, but it is along the same lines of phenomcnological models indicating that a less 
asymmetric DA is preferable to describe the nucleon form factors. Departures from the SU(6) symmetric model are 
necessary in the description of the neutron form factors [75| and the large x behavior of the neutron structure functions 
in deep inelastic processes [73j . Such SU(6) symmetry breaking contributions would lead to an asymmetric DA. This 
will be studied in a broader framework in a forthcoming paper [8fj| . In any case, after evolution from the low-scale 
of the model to Q 2 = 1 GeV 2 the first and second DA moments, calculated within the present model and shown in 
Table HH already compare well with lattice QCD results [3l| . 

In contrast to the nucleon DAs that have been studied for a long time, only very recently attention to the nucleon- 
to-meson TDAs has been drawn, and the possibility of having some information from experiment has been suggested. 
Here, for the first time a model calculation has been presented for the eight leading twist ./V — * tt° TDAs. They 
receive contribution from the fluctuations of the nucleon in (pir°) and (A + 7r°) subsystems and can be expressed as the 
convolution of the baryon DAs with the probability amplitude to find the corresponding baryon-meson component in 
the nucleon. The relative contribution of these components depends on the momentum transferred between the initial 
nucleon and the final pion as well as on the different spin configurations of the intermediate baryon. In particular, 
the A plays a special role in the case of the tensor TDAs which involve configurations with helicity ±3/2, while the 
interplay of the nucleon and A contributions with helicity ±1/2 determines the different shape of the vector and 
axial-vector TDAs. 

Results have been shown under kinematic conditions reachable, e.g., at GSI-FAIR as proposed in Ref. [H|, but the 
model can easily and will be applied to other kinematics such as those proposed to study at Jlab [56| . 
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Appendix A: SPIN COMPONENTS OF THE BARYON LIGHT-CONE WAVE FUNCTIONS 

ive the explicit results for the spin-dependent component of the LCWFs of the proton and 



In this Appendix we g: 
A state. 

In the case of the proton, we have 
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In the case of A we can have the following spin configurations: 
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For the A in the spin | state: 
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For the A in the spin —I state, indicated with Jj.: 
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Throughout Eqs. f&l j) - f&48]) we used the following definitions: a.; = (m + x^Mq), N(ki) 
k^ = fc^ x H~ ikiy, and fc,^ = /c^ ^ ikiy. 
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